An intrinsic route to melt fracture in polymer extrusion: 
a weakly nonlinear subcritical instability of viscoelastic Poiseuille flow 
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As is well known, the extrusion rate of polymers from a cylindrical tube or slit (a "die") is 
in practice limited by the appearance of "melt fracture" instabilities which give rise to unwanted 
distortions or even fracture of the extrudate. We present the results of a weakly nonlinear analysis 
which gives evidence for an intrinsic generic route to melt fracture via a weakly nonlinear subcritical 
instability of viscoelastic Poiseuille flow. This instability and the onset of associated melt fracture 
phenomena appear at a fixed ratio of the elastic stresses to viscous stresses of the polymer solution. 

PACS numbers: 47.20.Ft, 47.50.+d, 83.60.Hc 



If a polymeric fiber is produced by extruding it from 
a so-called die, a cylindrical tube or planar slit, the sur- 
face often exhibits undulations or irregularities at higher 
flow rates — see Fig [|. For increasing flow rates the 
undulations become progressively stronger, so much so 
that it can eventually cause the extrudate to break - 
hence the name melt fracture. A detailed understanding 
of this phenomenon has remained elusive for already over 
30 years. 

Depending on the geometry and type of polymer, var- 
ious types of phenomena seem to occur [Q ||, ||. The 
short wavelength deformations of the interface often ref- 
ered to as "sharkskin" instability appear to originate at 
the outlet: the extrudate quasi-periodically sticks to the 
outlet, widens, snaps loose and narrows. In the "spurt- 
flow" regime, the extrudate shows intermittent bands of 
smooth and irregular surfaces; there is good evidence 
that this has to do with a stick-slip phenomenon at the 
wall of the die. In spite of the multitude of possibilities, 
there is every reason to believe that when these instabil- 
ities are absent, as in the experiments of Fig. |l], polymer 
flow still exhibits some elusive generic bulk flow instabil- 
ity: According to the engineering literature ||, a quali- 
tative change in the flow behavior appears to occur at a 
more or less constant ratio of the normal stress difference 
of the melt over the shear stress for almost any polymer. 
Isolating this intrinsic mechanism is the main purpose of 
this paper. Our contribution also opens up new avenues 
for a more general understanding of the rich variety of 
complex fluid flow instabilities, such as polymer turbu- 
lence |j) or the fact that complex fluids such as granular 
media are prone to instabilities (shear banding) in stan- 
dard shear configurations. 

The first linear stability analysis of the flow in the die 
of a viscoelastic fluid described by the so-called Oldroyd- 
B constitutive equations [Q, performed already some 25 
years ago ||, showed that the flow was stable. Since then 
it is has been generally accepted that Poiseuille flow of 
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FIG. 1: Illustration of the surface irregularities that occur 
in the extrusion of a polymer fluid from a capillary tube 
(the wider structure at the top) in the absence of shark-skin 
or spurt flow instabilities. For small speeds, the surface is 
smooth (a), but beyond some flow speed surface irregularities 
develop, (b) and (c). The wavelength is of the order of twice 
the diameter d of the capillary. After M . 



viscoelastic fluids is linearly stable. The absence of a 
clear linear instability explains the focus on other mech- 
anisms, like those mentioned above. In this paper we 
present the result of an analytical nonlinear amplitude 
expansion which show that viscoelastic Poiseuille flow 
exhibits a weakly nonlinear (or "subcritical") instability 
due to normal stress effects; this instability appears to 
make melt fracture phenomena unavoidable for polymer 
fluids with normal stress effects, even if care is taken to 
suppress instabilities associated with the wall or shape of 
the die. Recent experiments on polymers whose rheology 
is well captured by the Oldroyd-B model confirm our 
scenario as well as our quantitative predictions. 

Although polymer fluids can exhibit all kinds of com- 
plicated relaxation phenomena, the basic feature essen- 
tially all polymer fluids (and in fact most complex fluids) 
share is the occurrence of elastic stress effects: when the 
shear rate is sufficiently strong that the polymers be- 
come stretched by the flow gradients, the forces along 
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the normals of a little cubical fluid element are different 
in different directions, unlike what happens for a New- 
tonian fluid where the pressure is isotropic. Microscop- 
ically, this effect is due to the fact that in sufficiently 
strong shear flows, the polymers get stretched (like little 
rubber bands) and acquire a nonisotropic orientational 
distribution. Thus they "pull" differently on fluid ele- 
ments in the direction along the flow and along the shear 
gradient. The simplest model to capture this effect (and 
which hence has become the working horse of theoretical 
studies of viscoelastic polymer flow [Q, ||, is the so- 
called Oldroyd-B model; we study it in the limit of large 
polymer viscosity, in which case it is refered to as the 
Upper Convected Maxwell (UCM) model. In this regime 
it is defined by the following constitutive equation for the 
shear stress tensor r in terms of the velocity shear tensor 
Vd through 
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A [dr/dt + v-Vt- (Vv) j t-t- (Vw) 
= -r)(Vv + (Vw) t ), 



(1) 



This constitutive equation is characterized by one single 
relaxation time A. The first two terms between square 
brackets together constitute the total time derivative of 
a fluid element moving with the flow; the other two non- 
linear terms are required by frame-independence, e.g. the 
fact that a solid-body rotation of the fluid does not lead 
to elastic stress effects. For fluids given by this consti- 
tute equation, the velocity profile in Poiseuille flow is 
parabolic, just as in Newtonian fluids. Moreover, in a 
cylindrical tube with radius R and coordinates (0, r, z), 
the above-mentioned Weissenberg number (also called 
Deborah number) which is defined as 



Wi = 
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becomes, upon following the usual convention to denote 
the shear rate dv z np / dr\ wa n by 7, 



Wi = 2A 7 = iv max X/R. 



(3) 



Note that the Weissenberg number defined in (g) is in- 
deed the ratio of the so-called first normal stress differ- 
ence r rr — t zz at the wall over the shear stress at the 
wall. Newtonian fluids are isotropic and hence the stress 
difference is zero; for Oldroyd-B model fluids, the normal 
stress difference T rr — r zz ~ 7 2 Q while the shear stress 
is linear, t tz ~ 7; as a result Wi is simply linear in 7. 

The weakly nonlinear instability scenario that our cal- 
culations imply is illustrated in Fig. ||a. Along the verti- 
cal axis we plot the relative perturbation of the basic flow 
profile (e.g. the relative shear perturbation at the wall, 
or the relative shear stress perturbation at the wall) as 
a function of Wi. Perturbations of the flow whose ampli- 
tude is bigger than the value given by the thick dashed 
line are unstable; their amplitude grows to some non- 
linear value indicated by the full line. This line denotes 
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FIG. 2: Summary of our main result for the stability of vis- 
coelastic Poiseuille flow in a cylinder in the zero Reynolds 
number limit, (a) Semi-quantitative sketch of the bifurcation 
scenario. See text, (b) Our results for the threshold values 
of the amplitude beyond which the flow is unstable; these 
curves correspond to the dashed branch in (a). Dots indicate 
the critical value of the shear rate at the wall (normalized to 
the shear in the unperturbed case), squares the critical value 
of the relative shear stress perturbation, multiplied by 10. 



nontrivial flow behavior with a characteristic wavelength. 
Note that the dashed line never touches the horizontal 
axis, in agreement with the fact that Poiseuille flow is al- 
ways linearly stable || , although the threshold amplitude 
becomes very small for large Wi. Furthermore, note that 
the full line merges with the unstable (dashed) branch 
(this is called a saddle-node bifurcation point in techni- 
cal terms); below the corresponding value Wi c , Poiseuille 
flow is nonlinearly stable to perturbations of any am- 
plitude. In this paper, we focus on our results for the 
cylindrical die, since this is the most relevant case. For 
slits similar results are obtained (hJ. 

Fig. ||b shows our analytical results for the (dashed) 
branch that marks the threshold amplitude beyond which 
the flow us unstable. As described in more detail below, 
our analysis is based on an expansion to third order in 
the perturbation amplitude; beyond some critical value 
Wi c « 5, we find that the cubic terms in our expansion 
lead to an instability, and from this we are able to cal- 
culate the threshold value. Clearly, our results are fully 
consistent with the scenario of Fig. ^a: the flow is stable 
to arbitrarily small (linear) perturbations, but perturba- 
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tions as small as a few percent in the shear stress at the 
wall already render the flow unstable for values of Wi 
around 10. The critical value of Wi c we estimate from 
our calculations is also in agreement with the very recent 
precise experimental data on a model UCM fluid g [□], 
and with the range where long polymers have been re- 
ported to show a change in the flow behavior Q . 

Actually, a number of observations already made us 
believe that the scenario summarized in Fig.^| was a vi- 
able one before we embarked on our analysis: (i) For 
increasing Wi, the linear stability actually becomes arbi- 
trarily weak (the damping of the linear modes decreases 
as 1/Wi). (ii) In the zero Reynolds number limit, the 
flow of an Oldroyd-B fluid in a Taylor Couette cell (two 
concentric rotating cylinders) is linearly unstable above 
some well-defined value of the Deborah number, which is 
analogous to the Weissenberg number j|, [|. This lin- 
ear instability is due to the fact that "hoop stresses" 
generally make flow along curved streamlines unstable 
H 0. However, recent experimental investigations |ll| 
have clearly demonstrated that the instability is "subcrit- 
ical": as sketched qualitatively in Fig. ||a, the nonlinear 
flow branch corresponding to roll-type patterns in the 
Taylor-Couette cell has been shown to extend down 
to about half the critical value where the linear instabil- 
ity occurs (the point where the dotted line intersects the 
horizontal axis) . The subcritical character of the insta- 
bility has been argued not to depend on the curvature 
of the streamlines which causes the linear instability in 
the Taylor-Couette cell. Thus, it is reasonable to assume 
that both in viscoelastic Poiseuille flow and in viscoelas- 
tic Taylor-Couette flow there is a subcritical instability, 
and that the only essential difference is that in the first 
case the dashed branch never intersects the horizontal 
axis, while in the second case it does. (Hi) The transi- 
tion to turbulence for Poiseuille flow of Newtonian fluids 
is also subcritical 

Returning to our analysis, it is based on deriving 
the first nontrivial nonlinear term in an amplitude ex- 
pansion for the perturbation to the velocity and shear 
stress fields of Poiseuille flow of an UCM fluid. Con- 
ceptually, this method goes back to the classic work 
fL4| in which the subcritical nature of the transition 
to turbulence of ordinary fluids was established; since 
then, amplitude expansions have become a standard 
tool in the field of nonequilibrium pattern formation 
p5j . We write the constitutive equation (|l|) and the 
Navier-Stokes equation (which in the zero Reynolds num- 
ber limit is linear) in cylindrical coordinates and then 
study the evolution of the amplitude A of a perturba- 
tion A(5v r , Sv z , Srgg, Sr rr , Sr rz , ST zz )e lkz ~' luJt , where the 
vector is normalized such that the shear perturbation 
dSv z /dr at the wall equals A frl6| ]. To first nontrivial 
order we then obtain an equation of the form 




FIG. 3: Values of the critical shear stress amplitude r c (fc), 
the ratio of the perturbation in the shear stress r rz (normal- 
ized to the value r"" p of the unperturbed state) at the wall 
above which the perturbation is nonlinearly unstable, as a 
function of the wavenumber k of the perturbation in the z- 
direction. 



To linear order in A this equation simply reproduces the 
term cu(k) of the dispersion relation of a single mode 
gifcz-jwt. this term is already contained in the old analy- 
sis of H . In particular, since we know that every mode 
k is linearly stable, lmoj(k) < for all k. 

The crucial new feature of our analysis consists of cal- 
culating the coefficient C3 explicitly — although this is 
technically demanding, the analysis is standard and con- 
ceptually straightforward. It will therefore be discussed 
elsewhere fic}| . In particular the real part of C3 is of im- 
portance for determining whether or not the flow is non- 
linearly unstable: if Rec3 < 0, then the nonlinear terms 
increase the damping of the amplitude and the unper- 
turbed state is, within this approximation, also nonlin- 
early stable. On the other hand, if Re C3 > 0, then the 
nonlinear term promotes the growth of the amplitude, 
and in particular amplitudes 



\A\ > A c 



'lmoj(k) 
Rec3 



(5) 



dA/dt = -iuj(k)A + c 3 \A\ 2 A + 



(4) 



grow without bound. Hence, in this approximation A c 
constitutes the critical amplitude beyond which the flow 
is nonlinearly unstable. The results presented in Fig. ^fa 
for A c are obtained directly from our results for the co- 
efficient C3. In dimcnsionlcss units and with our nor- 
malization, A c immediately yields the relative shear rate 
perturbation at the wall, necessary to make the flow un- 
stable. With a numerical factor that we obtain from our 
analysis, A c can be converted into a value for the criti- 
cal relative shear stress ratio r c , the perturbation in the 
shear stress beyond which the flow is unstable. 

As we already discussed in connection with Fig ||b, for 
Wi > 5, we find that there is some range of wavenumbers 
k for which Rec3 > 0, and hence for which the Poiseuille 
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FIG. 4: Plot of the width of the fc-band where the correspond- 
ing modes render the basic Poiseuille flow profile unstable. 

profile is nonlinearly unstable. In Fig. ^ we show for three 
values of Wi the relative critical shear stress amplitude 
r c (k) beyond which the flow is unstable, as a function of 
the wavenumber k. The further one gets above the value 
Wi c where the instability sets in, the wider the band of k- 
values is where we find an instability; the values plotted 
in Fig. ||b for r c correspond to the minimum values of the 
curves in Fig. |^, but as this figure shows for larger values 
of Wi, the bands are very flat so that the strength of the 
instability does not appear to depend very sensitively on 
k within the unstable band. 

In Fig. |] we plot the band in which modes are non- 
linearly unstable beyond some threshold amplitude. It is 
important to realize that k in this figure is the wavenum- 
ber of the modulation of the Poiseuille flow inside the 
tube: Since the flow changes upon exiting the tube, this 
is not necessarily the same as the wavenumber of the in- 
duced roughness modulation outside the die, as in Fig. |l|c, 
but the two can be connected by equating the frequency 
of the modulations measured at a fixed position. Our 
numerical results show that the mode with kR w 3.75 
which according to Fig. || emerges at Wi = Wi c w 5 
moves with a velocity of about v « 1.9 w av where u av 
is the average flow velocity. If the nonlinearities do not 
change this speed to much, this suggests a frequency / at 
onset of about 3.75 v/(2ttR) s» 1.13w av /i?. Further above 
threshold our amplitude expansion to cubic order neither 
yields the speed of the finite-amplitude modulated stress 
pattern nor the selected wavenumber. However, it is rea- 
sonable to assume that in the saturated nonlinear regime, 
the pattern moves with speed v av in the tube. From the 
band of unstable wavenumbers shown in Fig. ^, we then 
conclude that the frequency of the extrudate modulations 
should then lie in the range 0.3v av /R < / < 0.72v av /R. 

The fact that the critical Weissenberg number Wi c w 5 
which we find is in good agreement with the transition 
value noted empirically is already a strong indication 



that the instability which we have identified lies at its ori- 
gin. Further independent evidence comes from a series 
of new experiments which we have performed on a range 
of PVA Borax polymer solutions, whose rheological be- 
havior is well described by the Oldroyd-B/UCM model 
(constant shear viscosity, normal stress proportional to 
the shear rate squared). Our scenario implies that the 
first transition to melt fracture should be hysteretic, i.e. 
occur at a higher value of Wi upon increasing the flow 
rate than upon decreasing the flow rate; this is indeed 
observed in the experiments. Likewise, the wavelength 
of the roughness modulations is in good agreement with 
the theoretical analysis. We refer to j?J for details. 

In conclusion, we shown that normal stress effects 
drive viscoelastic Poiseuille flow unstable. We speculate 
that many nontrivial flow patterns of complex fluids that 
emerge under simple shear have a similar origin. 
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